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Abstract 



The equilibration of hot and dense nuclear matter produced in the central 
cell of central Au+Au collisions at RHIC (y/s = 200A GeV) energies is stud- 
ied within a microscopic transport model. The pressure in the cell becomes 
isotropic at t ~ 5 fm/c after beginning of the collision. Within the next 15 
fm/c the expansion of matter in the cell proceeds almost isentropically with 
the entropy per baryon ratio S/A = 150, and the equation of state in the 
(P, e) plane has a very simple form, P = 0.15e. Comparison with the statis- 
tical model of an ideal hadron gas indicates that the time t ~ 20 fm/c may 
be too short to reach the fully equilibrated state. Particularly, the creation 
of long-lived resonance-rich matter in the cell decelerates the relaxation to 
chemical equilibrium. This resonance-abundant state can be detected exper- 
imentally after the thermal freeze-out of particles. 

PACS numbers: 25.75.-q, 24.10.Lx, 24.10.Pa, 64.30+t 
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I. INTRODUCTION 



The assumption that strongly interacting hadronic (or rather partonic) matter, produced 
in nucleus-nucleus collisions at high energy, can reach a state of local equilibrium (LE) [l|-|3| 
is one of the most important topics in relativistic heavy-ion program Q. The degree of equi- 
libration can be checked by fitting the measured particle yields and transverse momentum 
spectra to that of the thermal model in order to extract the conditions of the fireball at the 
chemical and thermal freeze-out (see, e.g., [|5|-|P2| and references therein). Here the equilib- 
rium particle abundances, which correspond to a certain temperature, T, baryon chemical 
potential, /xb, and strangeness chemical potential, /is, can be determined. However, the 
analysis is complicated, e.g., by the presence of collective flow of particles and the non- 
homogeneity of the baryon charge distribution in the reaction volume. The volume of the 
fireball is also taken as a free parameter in the thermal model. Various non-equilibrium 
microscopic transport-, string-, partonic-, etc. models have been applied to verify the ap- 
pearance of at least local equilibrium in the course of heavy-ion collisions at relativistic and 
ultrarelativistic energies |IS-E3]. To reduce the number of unknown parameters and simplify 



the analysis it has been suggested pO,pT| to examine the equilibrium conditions in the cen- 
tral cell of relativistic heavy-ion collisions, simulated by one of the microscopic models. For 
these purposes we employ the microscopic Ultra-relativistic Molecular Dynamics (UrQMD) 
transport model, which nicely describe the available experimental data on hadron-hadron 
and nucleus-nucleus collisions in a broad energy range 



The central cell of symmetric heavy-ion collisions is a convenient system because the 
velocity of its center of mass is essentially zero. The cell should be neither too small nor 



too large. Previous studies at energies from 10.74 GeV (AGS) to 1604 GeV (SPS) pTJJgT 



have shown that the cubic cell with volume V = 5x5x5 = 125 fm 3 is well suited for the 
analysis. The aim of the present paper is to study the relaxation of hot nuclear matter, 
simulated within the microscopic model, in the central cell in Au+Au interactions at RHIC 
(y/s = 2004 GeV) energy. 

The paper is organized as follows. The conditions of kinetic equilibrium in the cell are 
discussed in Sec. |TJ. Section |T| describes the statistical model (SM) of an ideal hadron 
gas employed for the comparison with the microscopic calculations. The relaxation of hot 
nuclear matter in the cell to thermal and chemical equilibrium is studied in Sec. [TV| Finally, 
conclusions are drawn in Sec. [V|. 



II. KINETIC EQUILIBRIUM 

First, the kinetic equilibrium has to be verified. The collective flow in the cell should be 
isotropic and small, so it cannot significantly distort the momentum distributions of particles. 
As shown in Fig. [I], the average longitudinal velocity in the central cell at RHIC energies 
reaches its maximum at ~ 1 fm/c after the Lorentz-contracted nuclei have completely passed 
through each other. The longitudinal flow rapidly drops and converges to the developing 
transverse flow. The longitudinal collective velocity of particles becomes smaller than Vfi ow = 
0.15 c at t = 5 — 6 fm/c at RHIC, while the similar process takes t = 8 — 10 fm/c at the 
AGS and SPS energies. This gives a small correction, only (mjv vj- low /2) f& 7 MeV for the 
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nucleoli spectra, and less than a MeV for pions, whereas the characteristic temperatures are 
around 140 MeV. 

Velocity distributions of hadron species become isotropic and nearly Maxwellian quite 
soon pOU21| . Isotropy of the velocity distributions results in pressure isotropy. The pres- 
sure in the longitudinal direction in the cell, calculated according to virial theorem ||25| , is 
compared in Fig. |](a) with the transverse pressure. The time of convergence of longitudinal 
pressure to the transverse one decreases from 10 fm/c to 5 fm/c with rising incident energy 
from AGS to RHIC, respectively. Thus, the kinetic equilibrium in the central cell in Au+Au 
collisions at RHIC energy is reached at t ~ 5 fm/c. 



III. STATISTICAL MODEL OF AN IDEAL HADRON GAS 

To verify that the matter in the cell is in thermal and chemical equilibrium one has to 
compare the snapshot of hadron yields and energy spectra in the cell with the equilibrium 
spectra of hadrons obtained in the statistical model (SM) of an ideal hadron gas containing 
essentially the same number of baryonic and mesonic degrees of freedom |p3|| . We recall 
the procedure briefly: the values of energy density, e, baryon density, p B , and strangeness 
density, ps, extracted from the microscopic calculations in the cell, are used as an input to 
the system of nonlinear equations of the SM 

^ iC = ^E£f M (^B,Ps), (1) 

i 

PB iC = ^E^-^ M ( T ^B,Ps), (2) 

i 

p fc = l_J2Si-Ni u (T,HB,H S ), (3) 

i 

containing the baryon charge, B i: and strangeness, S i} of the hadron species i. This 
procedure enables one to determine temperature, T, baryon chemical potential, pb, and 
strangeness chemical potential, p§. If the set of three parameters (T, p^,ps) is fixed, all 
macroscopic characteristics of the equilibrated system can be determined unambiguously. 
For instance, the particle yields, Nf u , total energy, Ef u , and pressure, P SM ; are calculated 
within the SM via the Gibbs distribution function 



fi{p,mi) oc exp ( — \/p 2 + mf/T + fii/Tj, (4) 

where p and are the momentum and the mass of the hadron species i, respectively. 
Namely, 

^ M = 7^fW ff(P,mi)d 3 P, (5) 



(2tt^) 3 „ 

E? U = J (P 2 + m 2 Y /2 f( P , m*) d*p, (6) 

pSM ^jS^Iw^ 2fipM (7) 
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with g.i being the degeneracy factor. The entropy density s = S/V can be determined either 
via the distribution function (|) 

s=-e r f t p > m *) t in K p > m ^ - ^ d ^ ^ 

i (2tt^) Jo 

or directly from Gibbs thermodynamical identity 

Ts = e + P SM - p B p B - p s Ps- (9) 

The chemical potential /ij of the ith hadron is determined by its baryon and strangeness 
charge, /ij = p B P>i + P&Si- Thermal and chemical equilibrium is assumed to occur in the cell 
when the spectra of hadrons in the microscopic cell calculations become close to the spectra 
predicted by the SM. 



IV. RELAXATION TO CHEMICAL AND THERMAL EQUILIBRIUM 

At the isotropic stage the total microscopic pressure is close to the grand canonical 
pressure, P SM ; as shown in Fig. @(a). Both pressures converge at about t — 5 fm/c, which 
is chosen as a starting point for the comparison with the SM. The input, {e, p B , p$}, and 
output, {T, /ib, /is}, parameters are listed in Table | together with the microscopic pressure, 
entropy density, and entropy per baryon in the cell. The final time of the calculation, defined 
from conventional freeze-out conditions e ~ 0.1 GeV/fm 3 or p to t ~ 0.5po ||21|| , corresponds 



to t = 20 — 21 fm/c. One may see that in spite of different initial conditions the freeze-out 
time in the central cell in heavy-ion collisions at RHIC energies is similar to corresponding 
freeze-out times at AGS |20| and SPS [21] energies. In accordance with general estimates 



the baryochemical potential at RHIC energies is small while the temperatures are well 
above the anticipated temperature for the QCD phase transition, T « 160 ± 10 MeV. The 
entropy per baryon in the cell varies slightly after the beginning of the kinetic equilibrium 
stage. A comparison with s/p B at lower energies is shown in Fig. 0(b). The hadron- 
string matter in the central cell seems to expand isentropically with s/p B = S/A ~ 12 
(AGS), 32 (SPS), and 150 (RHIC). Note that the results of the simulations at AGS and 
SPS energies are intriguingly close to the entropy per baryon values extracted from the 
thermal model fit to experimental data, namely, (S/A) AGS rs 14 and (S/A) SPS « 36 M. It 
is most interesting to compare the predicted value (s/pb) RHIC = 150 — 170 to the upcoming 
RHIC data. Together, isotropy in the pressure sector, isotropic and nearly Maxwellian 
velocity distributions of hadrons, and almost isentropic expansion of matter in the central 
cell strongly support the idea that a hydrodynamic regime is established after a certain time. 
In the following, we extract the equation of state (EOS) of hadron-resonance-string matter 
from the microscopic simulations at the quasi-equilibrium stage of nuclear collisions, which 
is part of this hydrodynamic picture. 

The evolution of the pressure with the energy density is depicted in Fig. ||(a) for AGS, 
SPS, and RHIC energies in the central cell, respectively. Several interesting facts can be 
gained from this figure. First, the pressure drops linearly with the decreasing energy density 
for all three energies. It means that the ratio P/e remains constant for the whole time 
interval of the (quasi)equilibrium stage. Thus the equation of state in the (P, e) plane has a 
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very simple form: P = 0.12 e at AGS, and P = 0.15 e at SPS and RHIC. We see that despite 
the significant difference in the center-of-mass energy of the nuclei, colliding at \/s = 17A 
GeV vs. y/s = 200 A GeV, the ratio P/e in the central cell is saturated already at SPS 
energies. 

Figure [|(b) presents the evolution of the EOS in the (T, ps) plane. The baryon chemical 
potential in the cell increases with decreasing temperature at AGS and SPS energies, but 
remains almost constant (and small) at RHIC. Also, all three curves can be well approxi- 
mated by a simple linear dependence. This fact deserves explanation. Substituting P = as 
in Eq. (Q) and omitting the term psps, because the net strangeness density in the cell is 
quite small, we get 

rpJl_ = g(l + a) _ 
Pb Pb 

If the ratio e/pB is constant, then the temperature T would depend linearly only on the 
baryon chemical potential /ib, and vice versa. However, the energy per baryon slightly 
decreases within the considered time interval for all three energies. At the final stage it 
drops to approximately 80-90% of its initial value. Therefore, from Fig. 0(b) it follows that 
the evolution of temperature as a function of e/pb and pb is mainly determined by the 
change of energy per baryon, and not the baryon chemical potential. 

Energy spectra, dN/(4:7rpEdE), of 7r's, iV's, A's, A's and K's in the cell at the stage of 
kinetic equilibrium are shown in Fig. [| at two different times, t = 5 fm/c and t = 10 fm/c at 
RHIC energies. At t — 5 fm/c the spectra of baryons seem to be in a reasonable agreement 
with those of the SM, while the slopes of the meson spectra are steeper compared to the 
predictions of the statistical model. This means, that the apparent temperatures of mesons, 
especially pions, are lower than the temperature given by the SM. Moreover, the energy 
spectra of pions can be decomposed on two components representing low and high energy 
pions. At if: = 5 fm/c the fit by two exponents to the pion spectrum yields the temperatures 
Tfcu « 114 MeV and T% igh « 158 MeV, while at t = 10 fm/c the results of the fit are Tf ow « 
100 MeV and Tfi igh « 150 MeV, respectively. Although the apparent temperature of pions 
from the high energy tail of the energy spectrum is closer to the temperature given by the 
SM, TfJf = 171 MeV, one may conclude that the thermal equilibrium in the central cell is 
not reached yet. The slopes of baryon spectra in the microscopic calculations at t — 10 fm/c 
are also steeper than the SM slope. Since the energy density e is the same in both models, the 
lower temperatures of hadronic spectra in the microscopic case indicate that the hot matter 
in the cell is not in chemical equilibrium. Therefore, one might expect that the fractions 
of mesons and resonances in the UrQMD cell are overpopulated. These extra-particles 
consume significant part of the total energy and effectively "cool" the hadron cocktail. The 
reheating proceeds via the absorption of mesons in the processes like 7T7t — > p, pp — > nn, 
or ttN — > A, AN —> NN, etc. Thus, our next step is to study the time evolution of 
hadron abundances in the cell. Note also that the difference between the temperatures of 
meson and baryon spectra can be explained by the fact that baryons experience many more 
elastic collisions and interactions with the resonance production per particle, which drive 
the system toward thermal equilibrium, than mesons. The mean number of interactions 
per, e.g., nucleon or delta increases from 5 interactions at t — 5 fm/c to 10 interactions at 
t = 10 fm/c. In contrast, pions (kaons), which are readily absorbed and produced by the 
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(10) 



resonance- string hadron matter, suffer only 0.5 (1) elastic collision per particle at t = 5 fm/c 
and 1.8 (2.5) ones at t = 10 fm/c. 

Figure [| depicts the yields of the main hadron species in the cell within the time inter- 
val 5 fm/c < t < 19 fm/c. It is interesting that microscopic spectra of pions, which are 
underestimated by the SM in the central cell at lower energies ||21|| , converge to the SM 
predictions at t ~ 15 fm/c. Also, the statistical model overestimates yields of nucleons, 
lambdas, and kaons, while the yields of both baryon and meson resonances are reproduced 
quite well. Since the baryon number and the strangeness are conserved in strong interac- 
tions, where is the rest of the hypercharge, Y = B + S, in the UrQMD calculations in the 
cell? As seen from Table [TT], where the partial densities of baryons Rb and antibaryons R B 
in the microscopic and macroscopic model are listed (the net values pe = -Rb — Rb are listed 
in Table |I|), the total yields of baryons and antibaryons in the SM are larger than those of 
the UrQMD. The hadron- resonance-string matter in the cell is not in chemical equilibrium; 
that is why the density of antibaryons is 2-3 times smaller than the equilibrium values. The 
results on strangeness densities are listed in Table |IT| for the strange baryons, mesons, and 
their antiparticles, respectively. Here the net Strang eness is ps = Rs -Rg -Rs' -Rj|^ - 
In line with the previous observation, the SM predicts significantly larger abundances of 
both strange baryons and strange antibaryons at 5 < t < 11 fm/c in the cell, while the 
densities of strange mesons are pretty similar. It is worth noting here that strangeness is 



underpredicted by the UrQMD if compared to experimental data at SPS energies, see [27 



At i ~ 12 fm/c the situation turns around: microscopic yields of strange baryons become 
closer to those of the SM, but the predictions for strange mesons diverge. Note also that up 
to t = 13 fm/c the net strangeness of baryons in microscopic calculations is larger compared 
to the net strangeness of baryons in the SM at the expense of the net mesonic strangeness. 
After t = 13 fm/c the fractions of the net strangeness deposited in baryonic and mesonic 
sectors in the microscopic calculations coincide with the SM results. 

Finally, the ratios of hadronic abundances are studied (see Fig. g). Here the results are 
presented separately for non-strange and strange baryons and mesons. In the baryon sector 
the resonances dominate over the strange and non-strange baryons until the end of the 
simulations. This can be taken as an indication of the creation of long-lived resonance-rich 
matter. Recall that the formation of resonance abundant matter plays an important role in 



the evolution of the system produced in nuclear collisions at lower bombarding energies p8 



It has been suggested that the significance of the resonances, mainly A's, for the system 
development should be diminished with rising center-of-mass energy of the collisions. This 
assumption is not confirmed in our simulations. The fraction of baryon resonances is almost 
70% of all baryons in the cell at RHIC at 5 < t < 19 fm/ c, while at SPS and AGS the number 
of baryon resonances decreases from 70% to 35%, and from 60% to 25%, respectively. The 
meson fractions of resonances shrink within the time interval 5 < t < 19 fm/c from 60% 
to 30% (RHIC), 50% to 20% (SPS), and 40% to 15% (AGS). But at RHIC energies the 
hot hadronic matter in the cell as well as in the whole volume of the reaction is meson 
dominated. The mesons, baryons, and antibaryons carry 90%, 7%, and 3% of the total 
number of particles in the RHIC cell at t > 10 fm/c (cf. 85%, 14.5%, 0.5% at SPS and 50%, 
50%, 0% at AGS). The microscopic ratios for mesons (Fig. |6], right panels) seem to be very 
close to the SM ratios. However, one has to keep in mind that the temperatures given by the 
SM fit are 40-50 MeV higher than the apparent temperatures of meson species. Since the 
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freeze-out occurs at t ~ 21 fm/c in the central cell at RHIC energies, the matter in the cell 
is frozen before reaching thermal and chemical equilibrium. This circumstance significantly 
complicates the extraction of the chemical and thermal freeze-out parameters by means of 
the standard thermal model fit. 

But can the formation of the resonance-abundant matter be traced experimentally? To 
answer this question the rapidity distributions of baryon resonances are plotted in Fig. [7] 
for central (b = 3 fm) Pb+Pb collisions at Ei a b = 1Q0A GeV and Au+Au collisions at 
■\/s = 200A GeV, respectively. (The UrQMD predictions for other global observables at 
RHIC energies can be found in PP|). Here only those resonances that decay into ground state 
hadrons, i.e. no final state interactions, have been accounted for. The rapidity distributions 
of baryon resonances at SPS energies have a characteristic Gaussian-like shape, while the 
dN/dy distributions of those at RHIC energies are nearly flat in the rapidity interval \y\ < 
3.5 resembling the Bjorken scaling picture of nuclear matter expansion at ultra-relativistic 
energies More than 80% of the baryon non-strange resonances at RHIC energies are 
still A's (1232). One can see that the density of directly reconstructible baryon resonances, 
especially A's and A's + S's, per unit rapidity at RHIC energies is quite high. Therefore, 
the resonance-rich hadron matter produced in central Au+Au collisions at y/s = 200A GeV 
can be detected. 



V. CONCLUSIONS 

In summary, the microscopic transport model UrQMD is applied to study the equation of 
state of the hot meson-dominated hadron-string matter produced in the central cell (V = 125 
fm 3 ) of Au+Au collisions at -/s = 200 A GeV. After the restoration of the isotropy of pressure 
gradients, the hadron spectra in the cell are compared with those of the statistical model of 
an ideal hadron gas. It is found that the expansion of matter in central collisions proceeds 
with constant entropy per baryon ratio in the central cell, S/A = s/pb — 150. Since the 
S/A ratios for the central cell in A+A collisions, calculated at AGS and SPS energies, are 
very close to the ratios extracted from the analysis of the experimental data, the expected 
value of the entropy per baryon ratio at RHIC lies within the range 150 < s/ps < 170. The 
microscopic pressure in the cell is also close to the SM pressure. It shows a linear dependence 
on the energy density in the cell, P = 0.15 e, which is similar to the P(e) dependence in the 



central cell at SPS energies plfl . The obtained result is in accord with the EOS P = 0.2 e, 
derived for an ideal gas of hadrons and hadron resonances J50]. The temperature T SM in the 
cell at RHIC energies is shown to be nearly independent of the baryon chemical potential 
p B . 

The further comparison of the energy spectra and yields of hadrons with the SM predic- 
tions shows that the full thermal and chemical equilibrium is not reached even at the late 
stage of the reaction. This means that the times t ~ 20 fm/c may be too short for the relax- 
ation process. Particularly, the deceleration of the relaxation to equilibrium is attributed to 
the creation of the long-lived resonance-abundant matter. For instance, pions which are fre- 
quently absorbed and produced in various inelastic processes, including formation and decay 
of resonances, experience on average only two elastic collisions per particle. This appears to 
be insufficient to reach thermal equilibrium in the system of strongly interacting particles. 
In turn, inelastic collisions try to restore chemical equilibrium in the cell. Non-equilibrium 



7 



densities of strange and non-strange baryons, like TV's, A's, E's, etc, and their antiparticles 
are still lower than the equilibrium values. Amazingly, the yields of resonances are in accord 
with the SM values from the very early times t ~ 5 fm/c. Even the abundances of pions are 
equalizing after t > 15 fm/c. However, this result should be taken with a grain of salt. The 
fitting temperature of the thermal model is higher than the inverse slope parameters of the 
energy spectra of particles in the cell, i.e., the temperature of the chemical freeze-out will 
be overestimated by the SM fit. 

According to microscopic calculations, resonance-rich matter survives until the thermal 
freeze-out when the contact between the hadrons is lost. It remains a challenging task 
to verify the formation of long-lived resonance-abundant matter in heavy-ion collisions at 
y/s = 200A GeV experimentally. 
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FIG. 1. Time evolution of the velocities of longitudinal (solid lines) and transverse 
(dashed lines) flows of baryons (left panels) and mesons (right panels) in an asymmetric cell 
< {x,y,z} < 2.5 fm of central heavy-ion collisions at = 10.7AGeV (upper row), 

Eiab = 160^4GeV (middle row), and \/s = 200 A GeV (bottom row), respectively. Arrows indi- 
cate the times t cr needed for Lorentz-contracted nuclei to pass through each other. 
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FIG. 2. (a) The longitudinal (solid lines) and the transverse (dashed lines) diagonal compo- 
nents of the pressure tensor P in the central cell of heavy-ion collisions at AGS, SPS, and RHIC 
energies compared to the SM results (dotted lines). 

(b) Time evolution of the entropy per baryon ratio, s/p B = S/A, in the central cell with V = 125 
fm 3 in heavy-ion collisions at AGS, SPS, and RHIC energies, respectively. Solid lines denote mi- 
croscopic calculations with the UrQMD model, dashed lines show the predictions of the statistical 
model. 
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FIG. 3. (a) The evolution of pressure P and baryon density e in the central cells of the 
heavy-ion collisions at AGS, SPS, and RHIC energies. 

(b) The same as (a) but for the (T, /iB)-plane. Solid symbols correspond to the stage of kinetic 
equilibrium, open symbols indicate the preequilibrium stage. The hatched area shows the expected 
region of the quark-hadron phase transition. 
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FIG. 4. Energy spectra of N (circles), A (squares), ir (triangles up), K + (open stars), K~ 
(triangles down) and A (solid stars) in the central cell of Au+Au collisions at RHIC at t = 5 fm/c 
(upper panel) and t = 10 fm/c (bottom panel). Dashed lines show the predictions of the SM. 
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FIG. 5. The yields of main hadron species in the central cell of Au+Au collisions at = 200^4 
GeV as a function of time as obtained in the model UrQMD (circles) together with the predictions 
of the SM (stars). 
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FIG. 6. Time evolution of the hadron to resonance ratio R = H/H res in the central cell 
of Au+Au collisions at RHIC shown separately for baryons (left panels), antibaryons (middle 
panels), and mesons (right panels), as well as for non-strange hadrons (upper row), strange hadrons 
(middle row), and total hadron yields (bottom row). Circles denote the UrQMD predictions, stars 
correspond to the SM results. 
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FIG. 7. The rapidity distributions of baryon resonances in Pb+Pb collisions at 

Eiab = 160 A GeV (left panel) and in Au+Au collisions at y/s = 200 A GeV (right panel). Collisions 
are calculated in the UrQMD with the impact parameter 6 = 3 fm. 
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TABLES 



TABLE I. The time evolution of the thermodynamic characteristics of hadronic matter in the 
central cell of volume V = 125 fm 3 in central Au+Au collisions at RHIC (^/s = 200A GeV) energy. 
The temperature T, baryon chemical potential strange chemical potential us, pressure P, 
entropy density s, and entropy density per baryon density s/ps, are extracted from the statistical 
model of an ideal hadron gas, using the microscopically evaluated energy density e ce11 , baryonic 
density Pg 11 , and strangeness density pg e11 as input. 
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TABLE II. The partial baryon densities of baryons and antibaryons, R B / B , given by the 
microscopic model and obtained by the SM fit to the URQMD data, within the time interval 
5<i<19fm/cin the central cell of Au+Au collisions at RHIC energies. 
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TABLE III. The partial strangeness densities, R l s , of baryons, mesons, and their antiparticles, 
given by the microscopic model and obtained by the SM fit to the URQMD data, within the time 
interval 5 < t < 19 fm/c in the central cell of Au+Au collisions at RHIC energies. 
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